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Abstract 

We study the effects of plasma viscosity on the dynamics of primordial mag- 
netic fields by simulating magnetohydrodynamics in the early universe by appro- 
priate non-linear cascade models. We find numerically that even in the presence 
of large kinetic viscosity, magnetic energy is transferred to large length scales. 
There are indications, however, that the inverse cascade stops at a given time 
which depends on the magnitude of viscosity. For realistic viscosities we do not 
find equipartition between magnetic and kinetic energies. 
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There are many ways that magnetic fields could have been generated at various 
stages in the early universe. In the theory of galactic magnetism such fields are often 
referred to as "cosmological" or "primordial", because they are created well before 
galaxies are formed. These fields may play the role of an initial condition for the 
galactic dynamo, a mechanism that would amplify magnetic fields and convert kinetic 
energy into magnetic [|I|. 

There are two major problems when invoking primordial magnetic fields as possible 
seed magnetic fields for the galactic dynamo. One problem is the very small length 
scale of such magnetic fields. The horizon scale at the time of the electroweak phase 
transition is just a few centimeters, corresponding to about 1 AU at the present time. 
This is nine orders of magnitude shorter than the radius of typical galaxies. However, 
this view is too simplistic, because nonlinear effects inherent in the magnetohydrody- 
namic (MHD) equations can lead to a redistribution of magnetic energy over different 
length scales. Unfortunately, MHD in 3+1 exceeds the possibilities of present day com- 
puters. Therefore one has to resort to models which simulate the MHD equations. In 
a previous paper 0, hereafter referred to as Paper 1, we developed a fully relativistic 
3+1 d version of the so-called cascade model appropriate for MHD and found that 
an inverse cascade is operative, whereby magnetic energy is continuously tranferred 
to larger length scales. The other problem is that around the time of recombination 
photon diffusion becomes very large and could smooth out all fluctuations This 
may then also destroy the magnetic fleld 0. The purpose of this paper is to show that 
this too is too simplistic a viewpoint, and that nonlinear effects most likely prevent 
this from happening. 

The basic equations have been presented and discussed in Paper 1. We started out 
from the fully relativistic MHD equations in expanding (flat) space and showed that 
all the terms arising from the expansion can be removed by using rescaled quantities 
and conformal time. 



where tn is the Hubble time and R{tH) is the expansion factor. Starting from random 
initial conditions, we obtained turbulent velocity and magnetic flelds, very much like 
those in ordinary (nonrelativistic) decaying hydromagnetic turbulence. In order to 
study the effects of this kind of turbulence we adopted a simple cascade model that 
captures the qualitative features related to turbulent energy spectra. Such models [Q 
have been rather successful in predicting even subtle corrections to Kolmogorov tur- 
bulence due to intermittency effects p. (For a recent review, see ref. [0.) Our model 
is more general in that it includes the effects of magnetic flelds. Unlike nonmagnetic 
turbulence, in the presence of magnetic flelds there is an inverse cascade of magnetic 
helicity, which leads to a transport of magnetic energy to larger and larger scales ||^. 
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Using our cascade model (also called shell model) we found that the integral scale 

lo= J {2Ti/k)EM{k)dk I j EM{k)dk, (2) 

where EM{k) is the magnetic energy spectrum, increases with the Hubble time ap- 
proximately like t^jf^- 

In Paper 1 we considered the case where the kinematic viscosity u was equal to the 
magnetic diffusivity 77 (inverse electrical conductivity) which, in turn, was assumed to 
be small (correponding to a large magnetic Reynolds number). However, around the 
time of recombination the photon mean free path became very large and photon 
diffusion became very efficient in smoothing out virtually all inhomogeneities of the 
photon-baryon plasma 0. This process if often referred to as Silk damping, which 
corresponds to a kinematic viscosity z/ ~ (in natural units). We have computed 
numerically the evolution of magnetic and kinetic energy spectra in two completely 
different cascade models and studied the effects of viscosity. Our results are presented 
in Figs. 1 and 2 and the main point can be summarized as follows: in the cascade 
models magnetic energy is transferred to large length scales even in the presence of 
large viscosity. It seems likely that the same is true also in full 3+1 MHD. 

Let us first use the cascade model of Paper 1 to investigate the effect of very 
large values of u on the magnetic field. In this model velocity and magnetic fields 
are described by the variables f„ and 6„, representing the "collective" behavior at 
wavenumber In the original model, kn = r" with r = 2. Below we shall also 
consider a continuous version of this model with r — * 1, which was orginally studied 
by Parisi [|10| in a hydrodynamical context. We reiterate here the salient features of 
the model. We present the equations of motion in a slightly different form with the 
negative and positive hehcity states split, to facilitate comparision with the subsequent 
continuous model: 
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with similar equations with t>+ v~ and 6.^ ^ b~. The model describes an expanding 
radiation dominated magnetic universe {p = p/3). In Eqs. (^-(^ t is the conformal 
time, and the units are such that p (for details, see p). Eqs. (H)-(§) conserve 



energy and helicity in the absence of magnetic diffusion rj and plasma viscosity u, and 
+ and - are related to magnetic helicity. Thus, if u = rj = 0, the quantities 

E = EilKl' + \b-\' + + K\'), H = - \bn\')/kn, (5) 

are conserved. This is a direct consequence of Eqs. In the case of ordinary 

MHD b is just the magnetic field. In a cosmological setting with a flat universe, b is 
the magnetic field multiplied by R{ty, and t is conformal time. The conservation of 
helicity exhibited above is nontrivial in three, but not in two dimensions, so that in 
this sense the cascade model is three dimensional. 

There does not exist any proof that the cascade model and the standard MHD 
equations are equivalent. There are features which are similar in both cases: the 
equations couple many different scales (making it hard to predict anything a priori, 
especially when the Reynolds number is large and nonlinear effects are important), 
they have similar conservation laws, and the equations of motion are similar. Also, 
in the pure hydrodynamic case the cascade model equations have been compared to 
experiments and good results on intermittency have been obtained. In any case, the 
properties of turbulence have not been derived from first principles, and the cascade 
model is therefore an interesting (toy?) tool. 

For large diffusion coefficients the equations become stiff and it is therefore essential 
to solve for the diffusion term exactly. Using the identity 

^ + uklv^ = e-^'^"*! (vr.e-^'''-') (6) 
we solve equations of the form 

^ + uklv^ = iV„(t) (7) 
by a modified second order Adams-Bashforth scheme 

vn{t + St) = Kn{vn{t) + |<5t[3Ar„(t) - Kr,N{t - St)]} (8) 

where 

Kn = e-"^""". (9) 

We start from an initial condition that yields a magnetic energy spectrum similar 
to that found at later stages. Fig. 2 of Paper 1 suggests that the spectrum has 
developed an inertial range which is approximately constant, EM{k) = const, for 
kd < k < ko, where ko = 2n/lo is the wavenumber corresponding to the integral scale 
and kd ^ /cq the wavenumber of the diffusive cutoff scale. We adopt 30 wavenumber 
shells (1 < < with A^ = 30) and place the cutoff wavenumber at n = 27 and use 
7] = 10"^^. We put Vn = initially and compare the results for two different values of 
1/ (10-2 and 10^). 
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Figure 1: Magnetic and plasma kinetic energy spectrum as a function of the wave 
number k in the cascade model for small {u — 10~^) plasma viscosity (left) and large 
{u — 10^) plasma viscosity (right). The highest time is t — 10^, corresponding to a 
Hubble time 10^^ 

The resulting spectra are displayed in Fig. 1. The results show that Vn is generated 
from bn, but it gets weaker if u is increased. At small u and at large times there is an 
approximate equipartition of the magnetic and kinetic energies with Vn = bn = BnR'^, 
where Bn is the unsealed magnetic field. Since we are using units in which the initial 
effective energy density p + p = 4:p/3 = 1 it follows that f„ is the Alfven velocity 
Bn/ \J^f> of the equilibrium plasma. For large u equipartition is lost, which signals 
the breakdown of the pcrturbative approach. As we shall argue later, the large v case 
is appropriate for the very early universe. More importantly. Fig. 1 shows that the 
evolution of 6„ is not affected by f„ and z/, but is rather governed by ohmic decay, 
bn ~ exp(— r^/c^t). It may also be seen that in the case of large kinetic viscosity i/=100 
the velocity modes decay approximately according to t'„ ~ cxp(— z/A;^t). 

We shall now study the effect of a large u on b. We have done some more detailed 
calculations in order to investigate what is behind the behaviour shown in Fig. 1. 
The results are shown in Fig. 2, which shows that for a sufficiently large viscosity, 
the inverse cascade stops. With a magnetic enegy of 10~^, the cascade goes over five 
decades in the most unfavourable case i/=100. The inverse cascade stops at some low 
value of k, and the field just disappears by Ohmic decay. One can then ask, at which 
temperature does the inverse cascade stop? To estimate this, one studies the ratio 
between the non-linear terms and the Silk viscosity term, i.e. the Reynolds number. 
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Figure 2: The behaviour of the cascade model for u = 10^ (left) and u = 10^^ (right) 
at large times (the last time slice is t = 10^^). The figure shows what happens to an 
initially flat spectrum (which does not extend to k=0). Here rj = 10^^^. 

We start at the annihilation t/^ ~ 1 sec, where we assume that there are no velocity 
fluctuations. This is not true if the primordial field originates at the electroweak phase 
transition T ~ 100 GeV, since there will always be velocity fluctuations associated with 
the magnetic field. However, to allow for the optimal conditions for Silk diffusion, we 
assume that there are no velocity fluctuations at the annihilation. We then have to 
estimate the resulting Reynolds number. 
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where the viscosity vsuk/ R is effectively given by the photon diffusion length, 

1 0^ 

y ^ GeV^. 
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As mentioned above, the inverse cascade spans over five decades in k. Thus, the typical 
k which enters in Re in eq. ([To|) is given by 



k = 10~Vxt^™ 



(12) 



where xt'^"' is the spatial scale scale of the magnetic field at the annihilation. The 
fraction x is thus determined by the mechanism for producing the primordial field. 
Before annihilation, we assume that the viscosity is so unimportant that it can be 



ignored. This is why we start at annihilation, where Silk diffusion becomes very large. 
The inverse cascade stops when Re is of order one, i.e. when the diffusion equals the 
non-linear terms. We may therefore conclude that in the particular case displayed in 
Fig. 2 (z/=100), with v = 10~^, the cascade stops close to recombination. 

To verify that the above results are generic and not just a feature of the model 
adopted, we will now consider another model of the full 3+1 MHD [|T0|. Formally it 
may be obtained from ®-(^ by passing to the limit r — >1, but it may also be viewed 
as a completely independent model of MHD. Writing r = 1 + e we get 
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In addition to eqs. (|13]) and (|1^ there are two equations more, obtained by making 
the replacements 6+ ^ b~ and ^ v~ . These equations conserve the quantities 

E = J f (Kr + \v-\' + \b^\' + r\'), H = J ^(|6+r - \b-\'), (15) 

for z/ = =0, provided the following boundary conditions are satisfied for k — >0 and 
/c — i> oo, 

kv~{v^y 0, kv^b^b~ 0, ^^"6+ ^ 0, 

kv-^{v-y 0,kv-b+b- ^0,kv+{b-y ^0, (16) 



and 



v-{b+)^ 0, v+{b-y 0, b+b-v+ 0, b+b-y- -> 0, (17) 



respectively. If the conditions at infinity are not satisfied, this corresponds to "diffusion 
at infinity". 

We shall now simplify the continuous cascade model by using a scaling first intro- 
duced by Parisi [|1^ for the case of pure hydrodynamics without helicity conservation. 
If one considers the original discrete equations (§)-(§) in the absence of viscosity, and 
assume that we start at time t=0 with a primordial spectrum 



b+{k, 0) = b-{k, 0) = F, v+{k, 0) = v-{k, 0) = 0, (18) 



then it is easy to see that the discrete cascade model equations have a solution with 
the scaling form 



b+ik,t) = ¥'B{k^+H),h-{k,t) = ¥'M{k^+n), 

v^{k,t) = k^u{k^+H),v-{k,t) = k^m{k'+H). (19) 

These equations imply that Es/k"^^'^ is a function of k^^^t only. 

If we insert the scaling (|T^ in the continuous cascade model Eqs. (|T3p-(|T^), then 
we see that powers of k cancel out neatly, and we are left with the coupled equations]^ 

— u'{x) = 3{2p + l)u'm + 4(1 + p)xmu' + 2{1 + p)xum' — m'^ 

- 3(1 + 2p)MB - 4(1 + p)xMB' - 2(1 + p)xBM' + M\ (20) 

and 

- B'{x) = 3pmB - 3puM + 2(1 + p)xmB' + (1 + p)xBm' - (1 + p)xuM' 

- 2{l + p)xMu' + {I + p)xmM' - {I + p)xMm' , (21) 

as well as two equations where the interchanges B ^ M and u ^ m have been made. 
Thus, the original set of 2N (with n < 2N) coupled differential equations have been 
replaced by only four. In these equations 



X 



k'^n (22) 



is the scaling variable. 

It should be emphasized that the continuous cascade model is a priori "as good 
as" the discrete one, since in both cases the conservation equations are satisfied when 
there is no diffusion (except for possible diffusion at infinity). It may of course turn 
out that phenomenologically one type is better than the other. 

In Paper 1 we found that if one starts with a primordial spectrum p =3/2, then 
there is an inverse cascade, transferring energy from large to small k. The scaling 



eqs. (|l^) satisfy this: For a fixed value of the variable x in (^21) , the functions 
B, M, u, and m have some definite values. Thus, as time is increased, these val- 
ues remain the same if k is diminished in such a way that x remains constant. Thus, 
the scaling relations (|19D predicts an inverse cascade. 

This is also seen if we consider the integral scale (the "correlation length"), given 
by Eq. (|^), which measures the large scale structures in comoving coordinates. We 
may consider Iq as the expectation value of 1/fc. From scaling we therefore expect that 

/o ^tfd+p), (23) 



^We have here selected the phase of b and v to be -7r/2. With a dynamical phase the number of 
equations double from four to eight, since each b and v have a real and imaginary part. 
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Figure 3: The energy spectra in the continuous cascade model for (a) small (z/ = 10^^) 
plasma viscosity; in this case there is equipartition and magnetic and kinetic energies 
are equal; (b) large {y = 10^) plasma viscosity; absence of equipartition is apparent. 
In both figures the curves are, from right to left, for conformal times t=l,10,100 and 
1000, and r] = lO'^. 



where is the initial primordial spectrum. From this we see that for p=3/2 one 
obtains 

k ^ t'i. (24) 

where the power is in essential agreement with the previously found value p ~ 0.25. 

Since we do not commit ourselves to any specific model, let us consider an initial 
spectrum fc^. The integral scale then behaves as in (^31) in comoving coordinates, 
which means that in physical coordinates it goes like 
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(25) 



Thus, if there exists a model with p=0, in such a case the large scale structures are 
of the order the horizon. The case p=0 corresponds to a scale invariant primordial 
spectrum dk/k. 

The scaling relation is in general only valid in the "inertial range" where viscosity 
can be ignored, because the viscosity term is inconsistent with the scaling. There is, 
however, one exception, where scaling and viscosity are consistent, namely p=l. Here 
—u'{x) on the left hand side of eq. (pOD is replaced by 



u'{x) — uu{x). 



(26) 



Similarly, in eq. (|2lD —B'{x) is replaced by 



- B'{x) ~r]B{x). 



(27) 



10"' 



10"^ 



W 10"^ 



10"" 



/t=1101.88 



10-'^' 2 

10 



t=32487 



/' 1=1.03359 

t= 1.0000 



10^ 



10" 



10^ 



Figure 4: The scaled energy spectrum of the discrete cascade model at different times 
(neglecting viscosity) . 

This is the reason we have preferred to perform the numerical calculations in the 
scaled version of the continuous cascade model for the case p—1. We took as the 
initial conditions u = m = 0, M = B = 1 and as in the discrete case, considered two 
viscosities (z/ = 10~^ and u — 10^). The units are again such that initially p = po = 1. 
Magnetic diffusion was set to 77 = 10~^. The difference with respect to the discrete 
case is the shape of the initial spectrum of 6„. In the discrete case 6„ was initially 
given by whereas here we took 6„(A;,0) = k. The results are depicted in Fig. 3, 
where the inverse cascade can clearly be seen, virtually independent of viscosity. They 
also confirm that equipartition is lost at large u. 

We have also checked the scaling of b for the data presented in paper 1. This 
corresponds to p=3/2. The scaling variable is thus x = k^^'^t, where it is important 
that t is the conformal time; in terms of the Hubble time, the scaling variable is 
kt^jf. From its derivation, scaling is only valid in this case in the inertial range, where 
viscosity can be ignored. We see that within the fluctuations of the raw data, the 
predictions are approximately valid for large times (see Fig. 4). 

Our results very much suggest that in the real MHD, inverse cascade is operative 
and is essentially not affected by Silk damping, except very late and perhaps for very 
weak fields (which anyhow are not interesting). In our units, in the early universe 
i/ 1. Thus we may conclude that it is unlikely that there is equipartition in the very 
early universe. Our relativistic approach remains valid roughly until recombination, 
after which the plasma becomes matter dominated. (In the non-relativistic regime p ~ 
R~^, which effectively produces an extra term of the form —{R/R)v in the equations of 
motion. However, because photon diffusion is so large, this braking due to expansion is 
unimportant and conclusions presented above still hold.) Therefore, it seems plausible 



that small random magnetic domains of the very early universe may grow to large 
scale fields, irrespective of the Silk diffusion. 
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